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Abstract 



A general form for ladder operators is used to construct a method to solve 
bound-state Schrodinger equations. The characteristics of supersymmetry 
and shape invariance of the system are the start point of the approach. To 
show the elegance and the utility of the method we use it to obtain energy 
spectra and eigenfunctions for the one-dimensional harmonic oscillator and 
Morse potentials and for the radial harmonic oscillator and Coulomb poten- 
tials. 
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I. INTRODUCTION 



Algebraic methods applied to quantum mechanical problems have a long history [|TJ. 
The development of the supersymmetric quantum mechanics |2||| has contributed to en- 
large these approaches. Examples of this contribution is the concept of shape invariance 
introduced by Gendenshtein [|J], and the study of coherent states for one- dimensional po- 
tentials using a superalgebra formalism 0-17). 

Recently, Balantekin has shown that shape invariance has an underlying algebraic struc- 
ture ||. That approach was applied to studied barrier penetration and exactly solvable 
coupled-channels problems [|l(| . 

In this article we show that the algebraic structure presented by Balantekin can be applied 
as a generalized method to solve one- and three-dimensional Schrodinger equations exactly. 
When one introduces shape invariance added to supersymmetric formalism it is possible to 
obtain a general approach to construct generalized ladder operators. The one-dimensional 
harmonic oscillator and Morse potential, the radial equation for the three-dimensional har- 
monic oscillator and Coulomb potential are analyzed using this algebraic approach. Different 
forms of raising and lowering operators, not defined in the shape invariance context, have 
been proposed [0-TT3J for these potentials. 



The usual approach to solve Schrodinger equations using superalgebra formalism is by 



constructing a hierarchy of Hamiltonians [14]. For the method showed in this paper that 



construction is not necessary. The n-th eigenfunction for a given shape invariant potential 
can be deduced from the ground state wave function and the appropriated ladder operators. 

In the next section we introduce the formalism to construct the generalized ladder op- 
erators. In the section III the one-dimensional potentials are solved using this approach, 
as examples, the harmonic oscillator and the Morse potential are explicitly solved. In the 
section IV the radial Schrodinger equation is studied and the spectrum and eigenfunctions 
for the three-dimensional harmonic oscillator and Coulomb potential are obtained. 

II. THE GENERALIZED LADDER OPERATORS 



From a bound-state Hamiltonian as (h and 2m = 1) 

dx 2 



H = -^- + V(x) (2.1) 



with 

H* n {x) = £„f„(i) , (2.2) 

the supersymmetric quantum mechanics can be constructed if H can be factorized in the 
bosonic operators: 

A ± (x;a) = T-^- + W(x;a) , (2.3) 
dx 

where W(x;a) is well known as the superpotential, which is a function of the position 
variable and a set of parameters, a, that represent space-independent properties of the 
original potential V(x). If so, we obtain the first supersymmetric partner Hamiltonian 
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H + = H - E = A + A~ , (2.4) 

where E is the ground state energy eigenvalue. The second partner Hamiltonian is 

E_ = A'A + . (2.5) 

The Hamiltonians H + and H_ have the same energy spectra except the ground state of 
H + , for which there is no corresponding state in the spectra of H_. We can see that, from 
Eq. ( p.4|) the superpotential satisfies the Ricatti equation: 

W(x; a) - W'(x; a) = V(x) - E = V+(x) . (2.6) 

It was introduced in 1983 by Gendenshtein M the shape invariance condition 

R(a 1 ) = V-(x;a )-V+(x;a 1 ) , (2.7) 

where R(ai) is independent of any dynamical variable and a\ is a function of a^. Potentials 
which satisfy this condition are exactly solvable. However, shape invariance is not the most 
general integrability condition, once there are exactly solvable potentials that seen do not 



show the shape invariance property ||15|| . In this work we consider shape invariance involving 
translations of the parameters a: 

a! = a + r], (2.8) 

where 7] is the translation step. The other kind of relation between the parameters so far 
studied, scaling JIB] , is not broached here because only for the case of translations the 
potentials are explicitly known and can be written in a closed form in terms of simple 
functions, whereas for the case of scaling the potentials can only be written as a Taylor 
series [|[. In this way, it is possible to define translation operators T(ao) as 

T(a ) = exp (v-^j ( 2 -9a) 

T-^ao) = Tt(a ) = exp (-V^j ■ (2.9b) 

These operators only act on objects defined in the parameters space. 

In order to introduce the generalized creation and annihilation operators we compose 
the translation operators T and the bosonic operators A as 

B + (a ) = A + (a )T(a ) (2.10a) 
B4a )=T\a )A~(a ). (2.10b) 

These operators B present the necessary algebraic structure || to identify them as ladder 
operators. As such they are generalizations of the harmonic oscillator ladder operators. 

The method to solve Schrodinger equation with shape-invariant potentials is similar to 
a factorization method applied to the case of a harmonic oscillator potential. Thus, the 
ground state, ^ (x;a ), must obey 
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B^(a )^ (x; a ) = A (a )^ (x; a ) = , 



(2.11) 



or 



^o(x; Oo) oc exp (— J W(x;a )dx\ . (2.12) 



Eq. ( 2.12 ) is a relation between the ground state and the superpotential commonly found in 



the super symmetric approach PJT3| . The excited states are obtained by the repeated action 



of the creation operator on the ground state 

# n (x; oo) = [B+{a )) n M^ «o) • (2.13) 

At this point we want to stress that, since the operators T do not commute with any in- 
dependent object and the operators A do not commute with any x-dependent object then, 
the operators B act on objects defined in the dynamical variable space only by the bosonic 
operators, and on objects defined in the parameters space by the translation operators. 
The normalization constant that transform Eq. (|2.12|) into a relation of equality, in general, 



depends on the parameters a. Thus, this constant is also affected by the action of the 
operators B on the eigenfunctions. The normalization constant is not explicitly computed 
in this work, it can be included during the process of the use of the operators B Rj or 
determined from the normalization condition of the eigenfunctions. 

In order to determine the energy eigenvalues, it is important to understand the rule 
played by the remainder R(ai), Eq. ( |2.7| ), in this context. From the translation operators 
defined in Eqs. (|2.9| ), we can write 



R(a n ) = T(ao)R(a n ^)T\a ) , (2.14) 

where 

a n = a + nr] , (2-15) 



is a generalization of the relation given by Eq. ( |2.8j ). The superalgebra permits to obtain 
a simple relation between the energy spectrum and the remainders R(a n ). To see how it is 
possible we use Eq. (|2.14|) to write down 



R{a n )B + {a ) = B + {a )R{a n ^) . (2.16) 
From Eqs. ( 2.14|) and ( |2.16| ) we have the following commutation relation 

[H + ,(B + r}= (jTR(a k ))(B + r . (2.17) 



\fe=i 

Now we apply this equation on ^o(x; qq) to have 



H + {[B+(ao)] n Mx\ «o)} = E {[B + (a )] n *o(z; a o)} , (2.18) 



\k=l 



since ^o(x; a ) is the ground state of H + with a null eigenvalue. Thus, [B + (ao)} n ty (x; a ) is 
an eigenfunction of H + with eigenvalue e n = J2k=i R( a k)- Therefore we have for the energy 
spectrum 
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E n = Eq + e n = E + R(a k ) , (2.19) 

k=l 

where Eq can be obtained directly from the factorization of the original Hamiltonian, 
Eq. (£2D or Eq. @. 

We would like to close this section by emphasizing that this algebraic approach is self- 
consistent and it permits to determine the energy eigenvalues and eigenfunctions of a bound- 
state Schrodinger equation from supersymmetric and shape invariance properties of the 
system. In this way this approach is a method for exact resolution of that equation. 

III. ONE-DIMENSIONAL APPLICATIONS 

To illustrate the method discussed in the preceding section we first analyse one- 
dimensional shape-invariant potentials. 

A. Harmonic Oscillator 

This is the simplest example of a shape-invariant potential. In this sense, it is interesting 
in order to become the approach clearer. The original Hamiltonian is 

and the factorization is 

H + = A + A~ = ( — f- + ax] ( -f- + ax] = — ^ + a 2 x 2 - a , (3.2) 
\ ax J \dx J ax 2 

thus the superpotential is simply given by 

W(x;a)=ax. (3.3) 
The supersymmetric partner Hamiltonian is 

d 2 

H_ = A-A + = -— + a 2 x 2 + a, (3.4) 
dx 2 

and the shape invariance condition, Eq. (|2.7|), establishes that 

R(ai) = a\x 2 + a — a\x 2 + a\ . (3.5) 

Since R(a\) must not depend on the x-coordinate then ao = ai and a comparison of H + , 
Eq. ( |3.2|) , with the original Hamiltonian H, Eq. ( |3.1|) , results that, in this case, all parameters 
are equal to 1. So we get a n — ao — 1, r\ — 0, R{a^) = 2 and for the ground state energy, 
that is given by the constant term in Eq. ( |3.2|) , Eq — 1. 

The harmonic oscillator is the simplest case of a shape-invariant potential because the 
parameter step is null (rj = 0) and from Eqs. (|2.9| ) and ( |2.10| ) the generalized ladder oper- 



ators are simply the same as the bosonic ones (B± = A^). Thus, we recover the standard 
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approach of the ladder operators method applied to the harmonic oscillator. In this way, 
from Eqs. ( 2.11Q and Q2.19 ) we have 



B_(a )^ (x; a ) = A (a )^ (x; a ) = 



(x;a J oc e 



-x 2 /2 



(3.6) 



E n = E + J2 R(a k ) = l + £2 = 2n + l, (3-7) 

k=l k=l 



as it should be. 



B. Morse Potential 

Algebraic approaches, but not developed within the context of shape invariance, for one- 
dimensional Morse potential can be find in Refs. f|TT|,|T^,|r^] . In Ref. || it is explored the 
shape invariance property to determine the SU(1,1) structure of the algebra associated to 
that potential. 

Here we start with the original bound-state Hamiltonian for the Morse potential, 

Factorization of this Hamiltonian gives 

d 2 2 
= -— + A 2 (l-e~ x ) - 2aA + a 2 + A (2a - 1) e"* , (3.9) 

where the superpotential is 

W(x; a) = A (l - e~ x ) - a . (3.10) 
The supersymmetric partner of H + is 

H = A~A + = + A 2 (l - e~ x ) 2 - 2a\ + a 2 + A (2a + 1) e~ x . (3.11) 

The shape invariance condition, Eq. ( |2.7|) , for this potential can be written as 

i2(ai) = (2a - 2a x + 2) Xe~ x - 2 (a - a x ) A - a\ + a 2 . (3.12) 

In order to have R{a\) independent of the variable x, it is necessary that a\ = ao + 1 then, 
from Eq. (|2.8|) , rj = 1. A comparison between Eq. (|3.9| ) and Eq. (|3.8|) shows that ao = 1/2 
and Eq = A — 1/4. The shape invariance condition gives to R(a,k) the following form: 

R(a k ) = -2 (a fc _i - a k ) A - a\ + a\_ x = 2A - 2k . (3.13) 

For this case the energy spectrum, given by Eq. (|2.19|) , is 
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E n = E + Y, (2A - 2k) = 2A 



k=l 



n H — 
2 



!>2 
" + 2, 



(3.14) 



that is the exact result 12,13 



From Eqs. ( p. 10 ), the generalized ladder operators can be written as 



B + (a ) = A + (a )T(a ) = 
£L(a ) =T\a )A-(a ) 



dx 



+ A 1 - e~ x 



a 



cxp 



'_d_ 

da. 



^ ( 9ao) dx ^ ( 



a 



(3.15a) 
(3.15b) 



Substitution of Eq. ( glO|) into Eq. ( CT ), gives for the ground state eigenfunction the 
following expression 



^0(2^ a ) oc e a ° x exp 



A ( x + e 



(3.16) 



The excited states are constructed by the repeated action of the creation operator B + , 
Eq. ( p. 15a ), on the ground state Eq. (|3.16|) . For example, for the first excited state we have 



\E'i(x; ao) oc B + (a )^/o(x; a ) = exp [(—A + a + 1) x] exp 



-Ae" 



{2A - 2\e~ x - 2a - l}. 

(3.17) 



These same results for the Morse potential eigenfunctions [0 can be obtained by other 
exact methods, just remembering that do = 1/2. 



IV. THREE-DIMENSIONAL POTENTIALS 

The algebraic approach developed in the section II is not restricted to one-dimensional 
potentials. As a matter of fact, it can immediately be extended to radial equations for 
central force potentials in more dimensions. In particular, here we present the results for 
the radial harmonic oscillator and the Coulomb potentials in three dimensions. 



A. Radial Harmonic Oscillator 

The original Hamiltonian for a radial harmonic oscillator is 

d 2 2 1(1 + 1) . , 

Following the same development used in one-dimensional cases, the Hamiltonian given by 
Eq. (f4.1|) can be factorized, for each value of I, as 



H + = H-E t i ) = l--=L + r- — I l- + r 

\ dr r \dr r 



d l + l\ d l + l 
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where the superpotential is 



W(r; I) = r 



l + l 



The supersymmetric partner of H + , Eq 

d 2 



is 



A' A" 



dr 2 



+ r2+ (/ + l)(/ + 2) 



21 - 1 . 



For this case, the shape invariance condition, Eq. 

= t , + (jo+iKjo+2) _ 2in _ t _ 



, can be written as 
h(h + 1) 



(4.3) 



(4.4) 



(4.5) 



Note that here, the parameter a gives place to the angular momentum number I. Imposed 
that the right-hand side terms of Eq. ( J4.5| ) must not depend on the variable r, we obtain 
h = l + 1. Then, 77 = 1 and R{h) = 4, and from Eq. Q = 2/ + 3. As in the 
one-dimensional harmonic oscillator case, the function -R(Zfc) is constant. Thus, the energy 
spectrum for the 3-D harmonic oscillator is 



= Eg } + ]T R(l k ) = 21 + 3 + An . 



(4.6) 



k=l 



The eigenfunctions can be obtained from the generalized ladder operators: 



B + (l)=A + (l)T(l) 



d l + l 
_ _)_ r 



exp 
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i?_(0=T T (/)A-(/) = exp -- 



d 

— + r - 

dr 



\dlj 
l + l 



(4.7a) 
(4.7b) 



The minimum energy state for a fixed I {n = 0), is obtained by the substitution of Eq. (|4.7b|) 
into Eq. ( j2.11| ), giving 



(4.8) 



The other eigenfunctions (for n > 1) can be obtained by repeated action of the creation 
operator, Eq. (|4.7a ) on the function given by Eq. ( |4.8|) . As an example, we have for n — 1 



2r 2 -2I-3) r l+l e~ r2 l 2 . 



(4.9) 



These results are obtained by other methods, for instance, in Ref. |18 



B. Radial Coulomb Potential 



The radial Hamiltonian for the Coulomb potential is 

d 2 1 1(1 + 1) 

= 1 o 1 n , 



dr 2 r 



(4.10) 
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which is factorized as 



H+ = H — E, 



(0 



l + l 



d 



l + l 



dr 2(1 + 1) r I \dr 2(1 + 1) 



d 2 



1 1(1 + 1) 



where 



dr 2 r 



W(r;l) 



+ 



4(Z + l) 2 



Z + l 



2(1 + 1) r 

The supersymmetric partner of the Hamiltonian ( |4. 1 1| ) is 



H 



A~A A 



d 2 1 | (/ + l)(/ + 2) | 



1 



dr 



a(i + iy 



Here again the quantum number I is used as the parameter a (ao = Iq). Usinj 
invariance condition, we obtain 



(4.11) 
(4.12) 

(4.13) 
the shape 



l 



l 



4(/ 



4(/i + l) 2 



(4.14) 



In order to have R(h) independent on r we must have l\ — Iq + 1 and thus, rj 

1 1 



1 and 



R(h 



4(/ + l) 2 A(h + l) 2 
For the k-th parameter, the shape invariance condition gives 

R(h) 1 1 



(4.15) 



4(Z fc _! + l) 2 A(l k + iy 



(l + n) 2 (l + n + iy 



(4.16) 



The energy eigenvalue is obtained from the constant term in Eq. ( |4.11| ); 



1 



4(/ + l) 2 • 

Thus, the energy spectrum for this case is given by 



^ 1 



k=l 



1 



1 



4(n + / + 1) 



2 ' 



(l + n) 2 (l + n + 1) 2 
The generalized ladder operators for the radial Coulomb potential are 

l + l 



(4.17) 



(4.18) 



B + (l) = A + (l)T(l) - 
B_(l)=T\l)A-(l) 



d 1 

+ 



dr 2(1 + 1) 



( d s 
eW \-dl, 



d 



6XP (I) 
l + l 



dr 2(1 + 1) 



(4.19a) 
(4.19b) 



The lowest energy level eigenfunction for each / is obtained using Eqs. Q2.ll ) and (|4.19b ) or 
by the substitution of Eq. ( ELT2D into Eq. ( gig ) 



(4.20) 



The other eigenfunctions (for n > 1) are calculated by the action of the creation operator, 
Eq. ( |4.19a| ) on the above eigenfunction. As an example, for n = 1 we have the following 
eigenfunction: 



*?(r) oc5 + (I)^(r 



(0/ 



1 



2/-3|r m e - r W + V . (4.21) 



.1 + 2 l + l. 

These results for the radial Coulomb potential are also obtained by other exact methods 



T2|. Usually for this case, the principal quantum number adopted is N = n + I + 1. 



V. CONCLUDING REMARKS 

A general constructive method to solve bound-state Schrodinger equations for shape 
invariant potentials is presented. This algebraic method is based on a generalization of 
harmonic oscillator raising and lowering operators. The introduction of hierarchy of Hamil- 
tonians, usual in this kind of algebraic treatment, is no longer necessary for the present 
approach. Solutions for one-dimensional harmonic oscillator and Morse potentials and for 
three-dimensional radial harmonic oscillator and Coulomb potentials are explicitly obtained. 

Solutions for other shape invariant potentials can also be obtained by this approach. In 
particular, radial harmonic oscillator and Coulomb potential in arbitrary dimensions ]T9|J21]| | 
are simple extensions of the results appointed here. 

An interesting discussion is concerning coherent states for one-dimensional potentials 
and supersymmetric quantum mechanics [7[]. In this way one has a elegant and compact 
method to introduce ladder operators and define coherent states from these operators. 
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